Water Maser Emission from the Active Nucleus in M 51 
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ABSTRACT 

22 GHz water vapor 'kilomaser' emission is reported from the central region of 
the Whirlpool galaxy M51 (NGC5194). The red-shifted spectral features (Vlsr~ 
560 km s^^), flaring during most of the year 2000, originate from a spatially unre- 
solved maser spot of size < 30 mas {< 1.5 pc), displaced by ;< 250 mas from the nucleus. 
The data provide the first direct evidence for the association of an H2O kilomaser with 
an active galactic nucleus (AGN). In early 2001, blue-shifted maser emission (Vlsr~ 
435kms~^ ) was also detected. Red- and blue-shifted features bracket the systemic ve- 
locity asymmetrically. Within the standard model of a rotating Keplerian torus, this 
may either suggest the presence of a highly eccentric circumnuclear cloud or red- and 
blue-shifted 'high velocity' emission from a radially extended torus. Most consistent 
with the measured H2O position is, however, an association of the red-shifted H2O 
emission with the northern part of the bipolar radio jet. In this scenario, the (weaker) 
northern jet is receding while the blue-shifted H2O emission is associated with the ap- 
proaching southern jet. 

Subject headings: galaxies: nuclei — galaxies: individual (M51 (NGC 5194)) — masers 
— radio lines: galaxies 

1. Introduction 

The discovery of a thin molecular Keplerian torus surrounding a supermassive object at the 
center of NGC 4258 has proven that high-resolution imaging of water maser emission from active 
galaxies is a powerful method to investigate the dynamics of the circumnuclear region of AGN 
(Miyoshi et al. 1995; Herrnstein et al. 1999). About 20 water megamasers are known to date (e.g. 
Braatz, Wilson & Henkel 1997), all of which associate with Seyfert 2 or LINER nuclei. However, 
not all megamasers form part of a circumnuclear torus (e.g. Claussen et al. 1998; Gallimore et 
al. 2001; Hagiwara et al. 2001) and there are also galaxies that are characterized by much less 
luminous H2O emission. Towards the inner 40" of the starburst galaxy NGC 253 and the nearby 
spiral M51 (NGC 5194) weak (Ln^o ~ 0.1 - 1.0 Lq) water masers were observed, which were 
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classified as 'kilomasers' (Ho et al. 1987; Nakai & Kasuga 1988). Water masers with such inter- 
mediate luminosities are also observed in the Milky Way (e.g., W49N) and in the nearby galaxies 
IC 10, M 33, and M 82 and are associated with HII regions (Greenhill et al. 1990; Argon et al. 1994; 
Baudry & Brouillet 1996). 

While all H2O megamasers studied in detail were found to arise from the central few pc 
of active galaxies containing an AGN, it is an open question whether or not some of the H2O 
kilomasers can also be linked to the nuclear activity of galaxies. Are H2O megamasers only the 
most luminous members of a much larger but less conspicuous population of nuclear masers as 
suggested by Ho et al. (1987)? The H2O kilomaser source M51 is one of the nearest galaxies 
with an active Seyfert 2/LINER nucleus (Terashima et al. 1998) and is thus a suitable target. To 
test the hypothesis of Ho et al. (1987) we present single-dish monitoring observations with the 
Effelsberg 100-m telescope and imaging observations with the VLA. 



2. OBSERVATIONS 

2.1. MPIfR 100-m radio telescope 

The Effelsberg 100-m telescope, operated by the Max-Planck-Institut fiir Radioastronomie, 
was used on several occasions between June 1995 and January 2001. Between 1995 and 1997 
observations of the 22 GHz H2O line were made in a position switching mode with offsets of 300", 
employing a K-band maser receiver with Tgys ~ 75 K on an antenna temperature (T^) scale. The 
backend was an autocorrelator covering a bandwidth of 25 MHz with 1024 channels, yielding a 
velocity spacing of 0.33 km s~^ per channel. Since 2000 we used a dual channel K-band HEMT 
receiver with Tgys ~ 60 K (after averaging both channels) on a scale. The observing mode was 
dual beam switching with a beam throw of 121" and a switching frequency of ~lHz, providing 
baselines of high quality. The backend digital correlator consisted of eight modules (four for each 
orthogonally polarized channel), each of them providing 512 channels and a bandwidth of 40 MHz 
(1.05 km s~^ channel separation). Amplitude calibration was obtained by repeated measurements 
of 3C286 (Baars et al. 1977; Ott et al. 1994). The beam efficiency of the telescope was ~0.33 at 
22 GHz. Pointing measurements were made towards the nearby continuum sources DA 251 and 
3C280; pointing errors were <15". 
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2.2. Very Large Array 

The NRAO'^ Very Large Array (VLA) was used in its A configuration on January 23, 2001, to 
observe the 22 GHz H2O hne toward M51. The observations were made under good seeing condi- 
tions employing a single intermediate frequency band of width 12.5 MHz divided into 64 channels 
of width 2.63 km s"''' each. The band was centered at Vlsr= 560kms~^, which is ~90kms~^ 
rcdshiftcd to the systemic velocity of the galaxy. Considering filter roll-off at the edges of the 
passband, the usable velocity range was 495 - 625kms~^ . Given that at the time of our VLA 
observations H2O emission at velocities lower than the systemic velocity had not yet been detected 
and considering the constraints imposed by the VLA's spectral line correlator modes, our setup did 
not cover the VLSR-range of that low- velocity H2O emission. During the local sidereal time interval 
12^5 to 15^^ the nuclear position of M51 was observed in a sequence of 5 min. scans, interspersed 
with 2 min. scans of the calibrator source J1419-I-543 for which we used a position of (a, 5)j2000 
= 14'^19™46?5974, +54°23'14!'787, determined with sub-milharcsecond precision by Johnston et al. 
(1995). The resultant integration time on M51 is about 100 min. The flux density, S, scale was 
determined by observations of 3C 286 and is estimated to be accurate to within 10%, assuming that 
S{3C 286) = 2.54 Jy at 22.2 GHz. The data were calibrated and mapped in the standard way using 
the NRAO's AIPS. The synthesized beam has a FWHM size of 0{'102 x 0''096, corresponding to a 
linear scale of ~5pc at D '^lOMpc, and a position angle of — 7?6. Maps of the H2O emission were 
made for each velocity channel. These 'channel maps' have a typical la rms noise level of about 
4mJybeam~^. We also produced a map of the H2O emission integrated over the velocity range 
over which significant emission was found in the channel maps, i.e. between 538 and 590kms~^ . 
This map has an rms noise level of 2.8 mJy beam~^. 

3. Results 

Fig. 1 shows the spectra obtained between 1995 and 2001. The upper panel displays single dish 
profiles showing particularly strong maser emission in spring 1997 and 2000. The latter flare lasted 
until almost the end of the year. Most of the emission is seen red-shifted relative to the systemic 
velocity (V"sys,LSR = 469±5kms~^; Kuno et al. 1995). For the first time blue-shifted emission 
is also observed: Clearly seen at the end of January 2001, the strongest component at Vlsr = 
430.5 km s~^ is also visible in the Dec. 2000 spectrum. The blue-shifted emission covers a range of 
410 to 440 km s""^ or -60 to -30kms~^ relative to Vgys, while the red-shifted emission is observed 
between 538 and 592kms~^ or about -|-70 to -|~120kms~^ relative to V^yg. Integrated intensities 
determined for the velocity intervals 410 - 440 km s^-*^ and 538 - 592 km s"-"^ are not correlated. 
Although not shown, all seven epochs include the velocity range near 690kms~^ , where Nakai et 
al. (1995) detected emission in early 1995. No H2O feature is seen at their velocity in our spectra. 
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The lower panel shows the VLA spectrum taken on Jan 23, 2001. These are the first molecular 

data from M 51 taken with subarcsecond resolution. Both with respect to flux density and lineshape, 
the spectrum is consistent with the single-dish data from December 2000 and January 2001. The 
540, 566, and 578kms~^ components are all seen. Within the limits of noise, there is no missing 
flux. The emission arises from a spatially unresolved maser spot of size <, 100 mas at a position 
of (a,(5)j2ooo = 13^29'"52?7085 ± 0?0004, +47°11'42^'796 ± O'.'OIS. The quoted uncertainty is the 
quadratic sum of the error due to the thermal noise in the map and the upper bound of the 
systematic error caused by the fact that the calibrator source and M51 are at different positions 
in the sky. This upper bound, which dominates the error budget, was determined by producing a 
map of 3C286 by applying calibration solutions determined from J1419+543 scans and determining 
the offset between 3C286's centroid position in that map and its true position. Since the angular 
separation between J1419+543 and 3C286 is significantly greater than the separation between 
J1419+543 and M51, and since M51 was observed over a much wider hour angle range than 3C286, 
we consider the above offset as a conservative upper bound on the position error of the M51 H2O 
emission. The absolute position errors of J1419-I-543 and 3C286 are negligible (see e.g. Johnston 
et al. 1995; Ma ct al. 1998). 

In none of the channel maps (see Sect. 2.2) does the centroid position of the H2O emission 
significantly deviate from the position determined for the integrated emission. From the measured 
peak flux density and the beam size we obtain a lower limit to the maser's brightness temperature 
of 10^ K. Assuming Gaussian distributions for both the H2O emission and the beam, deconvolution 
yields a maser source size of <, 30 mas and a peak brightness temperature of >, 10^ K. 

To search for 22.2 GHz continuum emission, we made a map from a u,'u-database obtained 
by averaging the line-free channels. No continuum emission was detected at a la noise level of 

0. 52 mJybeam"^. 

4. Discussion 

4.1. H2O maser versus radio continuum peak position 

The most important result of our observations is the determination of an absolute position of 
the maser with a la accuracy of '~15mas. A comparison of this maser position with an 8.4 GHz 

VLA A-array continuum peak position ((a, (5)j2ooo = 13''29™52?7101 ± 0?0016, +47°11'42{'696 ± 
0^'026; Kaiser et al. 2001) shows that the water kilomaser arises from a hotspot --^100 mas (5pc at 
D = W Mpc) north of the radio continuum peak that marks the nuclear position of the galaxy (e.g. 
Grillmair et al. 1997, 1998). A comparison with a 5 GHz VLA A-array continuum map of slightly 
lower angular resolution by Crane Sz van der Hulst (1992) also shows good agreement; in this case 
the position of the continuum peak is (a, (^)j2000 = 13'^29°^52?706, -|-47°11'42''55 (no errors given), 

1. e. the maser is located ~25 mas cast and ~250 mas north of the compact continuum source. 
From even lower resolution 5 GHz VLA B-array and 15 GHz C-array data. Turner &; Ho (1994) 
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determined a continuum peak position of (a, 5)j2000 = 13''29"^52!711, +47°11'42''61; an absolute 
positional uncertainty of 100 - 500 mas is quoted. A comparison with our data suggests that the 
maser is arising from a position 25 mas to the west and 185 mas to the north of the continuum 
peak. Since the registration of one map to another is problematical when different calibrators 
and frequencies are involved, we should consider the quoted individual position errors with utmost 
caution. A better approach is to compare the three offsets derived: With an average offset of 
180ib80mas, this is first direct evidence for a kilomaser located in the vicinity of an AGN. 

4.2. Extended molecular complex, nuclear torus or nuclear jet? 

The nuclear region of M 51 is complex. CO J=l-0 and 2-1 intcrfcromctric data (Scoville et al. 
1998; Aalto ct al. 1999) show a ^2" sized centrally located molecular cloud with a bulk velocity 
that is red-shifted with respect to Kjys- Is this peculiar cloud responsible for the predominance of 
redshifted H2O maser features? This is not likely. The average velocities of the main CO sub- 
complexes are in the range 462kms~^ < Vlsr < 545kms~^ (Scoville et al. 1998) and do not cover 
the velocity range seen in H2O (Fig. 1). At Vi,sR ~ 565kms~^ , the characteristical H2O velocity of 
the red-shifted features, Scoville et al. (1998) find CO J=2-l emission l" south-west of the nucleus, 
well outside the positional error bars of our data (Sect. 4.1). According to Scoville et al. (1998; 
their Sect. 3) the asymmetry in the CO distribution is significant, with the bulk of the emission 
being red-shifted with respect to the systemic velocity and arising from western offsets with respect 
to the nuclear position. 

HCN J=l-0 line emission is weaker but traces gas of higher density and may thus be a more 
appropriate tracer of the warm dense molecular gas that gives rise to 22 GHz H2O emission. In the 
HCN J=l-0 line, an asymmetry favoring red-shifted over blue-shifted emission (as seen in CO) is 
not apparent. At the 3(t level, Kohno et al. (1996) find HCN J=l-0 emission at 551-582 km s~^ 
in the inner 5" of M51 and identify a dense nuclear torus with R ~70pc and an apparent T^ot 
= 16kms^^ at its outer edge; the estimated inclination is i ~ 50° - 80° (much higher than the 
inclination of the large scale disk). Red-shifted emission with respect to V^ys is seen to the west, 
blue-shifted emission arises east of the nucleus. With an elongation along P. A. ~ 70°, the rotation 
axis (P.A. ~ 160°) is aligned with the nuclear jet. If the nuclear torus shows Keplerian rotation and 
if the very center dominates the mass ('-^lO''' Mq), V^ot = V - Vsys ~ 95kms~^ (for the H2O line 
velocities, see Sect. 3) should be reached at i? 2pc (40 mas). In view of the positional error bars 
(Sect. 4.1) this is not inconsistent with the observational data, but we note that H2O emission from 
V — Vsys ~ 95kms~^ would then arise from a location 40 mas west, not (as our data more likely 
suggest) ^180 mas north of the nucleus. Assuming that the blue- and red-shifted masers originate 
from those parts of the torus that are viewed tangentially, the asymmetry in V - T^ys (Sect. 3) 
infers a dense warm torus that is radially extended or the existence of highly eccentric orbits. The 
position offset between red- and blue-shifted features should be detectable with the A-array of the 
VLA. 
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Instead of inferring that H2O is tracing a nuclear torus, we may postulate that the emission 

is associated with the radio jet. As already mentioned, the radio jet and the rotational axis of 
the inner (HCN) torus are characterized by P. A. 160 - 165°. The maser position determined in 
Sect. 4 is not far from this axis (for the line connecting the 8.4 GHz continuum peak (Kaiser et al. 
2001) with the maser position, P. A. = 174° it 10°). Agreement with measured positions is thus 
better than in the Keplerian disk scenario (for the 8.4 GHz peak, the distance to the jet axis is 
15±35 mas). In view of unifying schemes predicting a Doppler boosted approaching and a Doppler 
dimmed receding jet, it is likely that the southern jet is approaching, while the northern jet is 
receding. Within this context the red-shift of the H2O emission is readily explained: Entrainment 
of the molecular gas by the jet may cause the velocity offset from systemic. 

It is well known that 22 GHz H2O emission is often associated with shocks (e.g. Elitzur 
1995) that may raise not only kinetic temperatures but also H2O abundances. Such shocks can 
be triggered (on small scales) by jets arising from young stellar objects or (on larger scales) by 
nuclear radio jets. Megamasers not oriented perpendicular but along the nuclear jets have already 
been observed towards NGC1052 (Claussen et al. 1998), Mrk348 (NGC262; Peck et al. 2001), 
NGC 1068 (the 'jet masers'; Gallimore et al. 2001), and Circinus (Greenhill et al. 2001) and may 
either arise from a fortuitously placed foreground cloud or, more likely, from a shocked dense region 
near the boundary of the jet. According to this latter scenario, the blue-shifted features from M51 
(see Fig. 1) should be associated with the southern jet. This can be checked observationally, since 
position offsets should be large enough to be detected by VLA A-array observations. 

4.3. Conclusions 

Presenting the first direct evidence for a nuclear H2O kilomaser, our data are not inconsistent 
with the suggestion of Ho et al. (1987) that there must be a family of weak nuclear 'megamasers'. 
The kilomaser in M 51 shows characteristic properties of megamasers. Only the measured luminos- 
ity is smaller. Interpreting slight offsets between nuclear and maser position and postulating that 
the northern jet is receding, we conclude that the masers in M51 are most likely associated with 
the nuclear jet. Positional uncertainties are, however, too large to definitely rule out emission from 
a circumnuclear disk. 

We wish to thank W. A. Baan for providing a high precision J 2000.0 8.4 GHz continuum peak 
position of M51, W. A. Altenhoff for the selection of a sufficiently accurate program to convert 
coordinates from B 1950.0 to J 2000.0, and A. B. Peck for critically reading the manuscript. We 
also appreciate P. T. P. Ho for his suggestions on the manuscript. 
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Fig. 1. — Upper panel: 22 GHz H2O spectra taken with the Effelsberg 100-m telescope between 
June 1995 and January 2001 from M51. Note the enlargement of the temperature scale by factors 
of 2 or 3 for the most recent spectra marking the end of the flare. For the first two epochs channel 
spacings are 1.3kms~^ ; for the other spectra the spacing is 1.05 km s~^ . la noise levels are, in 
chronological order, 22, 17, 28, 15, 15, 13, and 9mJy. Lower panel: VLA spectrum from January 
23. For details, see Sect. 2.2. 
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